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ON OVERLOAD IN A STORAGE MODEL, WITH A 
SELF-SIMILAR AND INFINITELY DIVISIBLE INPUT 

By J. M. P. Albin 1 and Gennady Samorodnitsky 2 

Chalmers University of Technology and Cornell University 

Let {X(t)}t>o be a locally bounded and infinitely divisible stochas- 
tic process, with no Gaussian component, that is self-similar with 
index H > 0. Pick constants 7 > H and c > 0. Let v be the Levy mea- 
sure on R [0 ' oo) ofX, and suppose that R(u) = v{{y £ ft' ' 00 ) : sup t > y(t)/(l + ct 1 ) > 
is suitably "heavy tailed" as u — > 00 (e.g., subexponential with posi- 
tive decrease). For the "storage process" Y(t) = sup s>t (X(s) — X(t) — 

C (s ty), 

we show that P{sup ae r t ( u ji Y(s) > it} ~ P{Y(t(u)) > u} as it — > 00, 
when < £(w) < i(u) do not grow too fast with u [e.g., t(u) = o(u 1 >'~')]. 

1. Introduction. Let X = {X(t)}t>o be an infinitely divisible (i.d.) stochas- 
tic process, with no Gaussian component, that is self-similar with index H > 


(fl-s-s.). 

Given constants c> and 7 > H, we consider the storage process 
(1.1) K(t) =sup(X(s) -c(s-t) 7 ) fort>0. 

Intuitively, an H-s.s. process grows as t H with time t, and so 7 > -ff should 
make Y finite valued. Nevertheless, this is not so in general (see Example 
2). The assumptions in our theorems will, however, ensure such finiteness of 
Y. The reason for the name "storage process" comes from the case 7 = 1, 
with X(t) denoting the total inflow into a storage facility by time t, and c 
the (demand) rate at which stock at the facility is depleted; then Y(t) tells 
how much extra storage capacity one will need in the future over what is 
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being used at time t. For an input process X with stationary increments 
(si.), the storage process Y is stationary (if finite). 

The process Y has been used in financial applications under the name 
of "drawdown" [e.g., Dacorogna, Gengay, Miiller and Pictet (2001)], and is 
important in queueing applications; for example, to model teletraffic, when 
X is Gaussian if-s.s.s.i. (i.e., fractional Brownian motion), with H > ^ and 
linear service 7 = 1 [e.g., Norros (1994) and Piterbarg (2001)]. In this case, 
building on Hiisler and Piterbarg (1999), Piterbarg [(2001), Theorem 5] gave 
a version of the remarkable property (1.2), that was a triggering influence 
for us. Recognizing this, we name that property after him. However, the 
Gaussian problem Piterbarg studied is very different from ours with i.d. 
processes, and his proof, by Gaussian field theory, does not relate to non- 
Gaussian settings. 

In Section 3 we study the probability for overload during a time interval, 



P< sup Y(s) > u > asii^oo. 
Lefo.tl J 



sG[0,t] 

If for each choice of a constant t > 0, it holds that 

P{sup sg r t i Y(s) > u} 

(1.2a) lim — ' = 1 whenever < t(u) < t, 

v > u^oo P{Y(t(u))>u} ~ V J ~ 

then we say that the process Y has the Piterbarg property. The similar 
statement, for which the length t = t(u) of the interval may depend on the 
level ti, 

(1.2b) lim P { sup *g[ '5")] y ( s ) > = 1 whenever 0<t(u)<t(u), 
K > u^oo P{Y(t(u))>u} - \ J - \ j, 

will be referred to as the generalized Piterbarg property. 

One indication of the unusual behavior of Y is that (1.2) implies, as 

u — ► 00, 



sup Y(s) > u 
se[o,t] J 




P{Y(U) > u} 
P{sup 8e[0jt] Y(s) > u} 



1. 



Thus overload periods within [0, t] are long enough to include any t±,...,t n £ 
[0,*]. 

This last conclusion leads us naturally to the question whether one can 
replace the minimum /\™ =1 Y(ti) taken over a finite collection of points in 
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[0,t] by the infimum over the entire interval. That is, we would like to know 
if 

M . . ,. P{sup se[0 , t ] Y(s) > u} 
(1.2c) lim „ r . - — ' . r = 1. 

n-^oo P{mf se[0)t] Y(s) > u} 

This we call the strong Piterbarg property, whether or not t is a function of 
u. 

With v being the Levy measure on IRl ' 00 ) of X (see Section 2.3), denote 

(1.3) i?(u) = lAi/(yeR (0 '°° )nQ : sup > u\) for «6R. 

VI te(o,oo)nQ 1 + c * 7 )J 

We will make assumptions about "heavy tails" for the function R (e.g., 
subexponentiality together with positive decrease; see Section 2.1). Under 
additional technical assumptions on X, we establish the generalized Piter- 
barg property, when t(u) does not grow too fast with u [e.g., t(u) = o( , u 1 / 7 )]. 
Under the same assumptions we will show that the strong Piterbarg prop- 
erty holds as well. Under certain weaker assumptions, we prove a weaker 
so-called O-version of (1.2), that is, that the probability ratios in (1.2) are 
bounded away from zero and infinity. 

Our main "external tool" in proofs is Theorem 2.1 on subexponential func- 
tionals of i.d. processes by Rosihski and Samorodnitsky (1993); see Section 
2.3. 

The Piterbarg properties are quite unusual. For example, only a degener- 
ate 

a-stable or Gaussian process Y can have them; see Example 7. 

In Section 4, we give a discussion, with examples of application, and coun- 
terexamples, for i.d. H-s.s. processes X given as stochastic integrals with re- 
spect to heavy-tailed i.d. random measures (see Section 2.3). This includes 
a-stable processes. 

2. Classes of functions and stochastic processes. It will be convenient 
to devote a separate section to describe classes of functions and stochastic 
processes, that feature in the rest of the article. In addition, some basic rela- 
tions between these classes, and some important representation properties, 
are listed for easy reference. 

2.1. Classes of functions. In this section, /:R— > (0, oo) denotes a non- 
increasing function with lim u ^. QO /(u) = 0, and g:M — > (0,oo) a measurable 
function. 

The function / is subexponential, f £ S, if there exist independent iden- 
tically distributed random variables £ and r], such that 

/(u)~P{£>u} and P{£ + rj > u} ~ 2P{£ > u} asu^oo. 



4 



J. M. P. ALBIN AND G. SAMORODNITSKY 




lim sup < 1 for some A > 1 . 

u^oo flu) 





Note. The definitions of OR, PD and ER are more complicated than 
those given above for a general nonmonotone /. 

The function g is regularly varying with index p G R, g G RV(p), if 

lim ^ ) = \ p for A > 1 (or, equivalently, for A > 0). 

Here the convergence must, in fact, be locally uniform. 

Notice that the function go log belongs to RV(0), which we denote g £ C, 



lim — — — = 1 for A > (or, equivalently, for A G R). 

We have ER n PD C ORnPDnS, and a monotone / G U P<0 RV (p) 
belongs to all these classes. Further, OR n C C 5 C £. 

The classes of functions above, and the listed relations between them, 
are well known from the literature. See, for example, Bingham, Goldie and 
Teugels (1987). 

2.2. Classes of stochastic processes. In the remainder of this article, 
X = {X(t)}t>o denotes a separable stochastic process, that is continuous 
in probability and locally bounded (bounded on any given compact inter- 
val) a.s., and is denned on a complete probability space P). We refer 
to these requirements as Condition X. Depending on the context, further 
requirements on X will be imposed later. 

We write {X(t)}t>o = X when the finite-dimensional distributions (f.d.d.'s) 
of the processes X and X agree. For example, X is stationary, if X(- + 

h)=Xfarh> 0. 

The process X has stationary increments, if 



if 



g(u + A) 



X(- + h)- X(h) = X - X(0) for each h>0. 
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The process X is self-similar with index H > (H-s.s.), if 

a~ H X(a-) = X for each a > 0. 

The process X is infinitely divisible (i.d.), if for each n G N, there exist 
independent processes {Xi(f)}t>o, • • • > {A^,(£)}i>o> such that 

X x = ■■■ =X n and Xi + • • • + X n = X. 

The process X is a-stable, a G (0,2], if for each n G N, there exists a 
constant process C n , such that, taking independent copies {Xk}^ =1 of X, 

n- 1 / a {X l + ---+X n ) + C n = X. 

In particular, it turns out, the process X is Gaussian if and only if it is 
two-stable. 

The process X is strictly a-stable if, taking independent copies {-Xfc}j£Li 
of X, 

n- 1/a (X 1 + ■■■ + X n ) = X fornGN. 

A process X is H-s.s. if and only if the Lamperti transformed process 
e~ X(e') is stationary [Lamperti (1962)]. 

a-stable processes are i.d. Clearly, an a-stable process X is strictly a- 

stable if it is symmetric a-stable (SaS) (a-stable with X = — X). 

Of course, the classes of processes mentioned above are all quite basic, as 
are the indicated relations between. See, for example, Samorodnitsky and 
Taqqu (1994) for further information, and for an extensive bibliography. 

2.3. I.d. stochastic processes. The f.d.d.'s of an i.d. process {X(t)} t £T, 
T = [0, oo), with no Gaussian component, can be described by means of a 
Levy measure v on the cylindrical <r-algebra B on R T , and a localization 
parameter \x e M T . 

Let 7r r be the projection of M T on W , and let B T be the Borel sets in 1R T , 
for tgT={tCT:1< #r < oo}. According to Maruyama (1970), a mea- 
sure v on B is a Levy measure for X, if v o 7T" 1 is a Levy measure onB T [i.e, 
if 1 A | • | 2 G L 1 (R T , v o 7T" 1 )] for each r G T, and there exists a // G M T such 
that 

E{e^' x >} = exp(i(0, fi) + [ (e i{e ' x) - 1 - i(6, k(x))) du(x)\ 
(2.1) ^ ^ RT ' 

for 6»GE (T) . 

Here we use the notation 

R( T ) ={xG]R T :#{tGT:x(t)/0} < oo}, 

( x , y) = ^ x(t)y{t) for x G M (T) and y G M T , 

= a?(*)l[_i i i](|a?(*)|) for x G E T and t G T. 
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[A general i.d. X can be represented as X = X\ + X 2 , with X\ and X2 
independent, X\ i.d. with no Gaussian component as in (2.1), and X 2 zero- 
mean Gaussian 

E{ e *<W} = exp(-± J2 e(s)em{x 2 ( s )x 2 (t)}\ 

for fle#>.] 

We now turn to the task of constructing and representing i.d. processes. 

Let (S, S, A) be a c-finite measure space, and put 60 = {^4 £ 6 : A(A) < 
00}. An (independently scattered) i.d. random measure (with no Gaussian 
component), with control measure A, is a map M : &q — > L°(f2, J) such that, 
for Ae&o, 

E { e ^M(A)| = exp f J ( i9m+ f ( e ite _ ! _ p (., dx)) dA) 

(2.2) Ua\ Jr J J 

for 0£R. 

Here the localization m G L°(5) satisfies {lyi^lyieeo L 1 (S, A), while p(s, •) 
is a Levy measure on M for s£S, such that p{-,B) G L°(5) for Borel sets 
SCR, and 

(2.3) F(A x •) = / p(s, -)A(<is) is a Levy measure on R for each ^4 G ©o- 

/a 

The stochastic integral J s f dM is well defined (in a P-sense) , for / G 
L°(5) with 



(2.4) 



(1 A \xf\ 2 )p(-,dx) d\ 



V 



mf+ (n(xf) - fn{x))p(-,dx) 



dX < 00 



E 

(2.5) 



[Rajput and Rosinski (1989), Section 2]. In that case, J s fdM is i.d., with 

h( ie L fdM )} 

= exp^J^iemf + Jj,e i6xf - 1 - ieK(x)f)p(-,dx)^J dx\. 

In the language of (2.1), the (process consisting of a) single i.d. random 
variable J s f dM has Levy measure v on ]R, and localization parameter p G M, 
given by 

u(B) = F({{s,x) G 5 X R:xf(s) G B}), 
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In particular, for example, by Feller [(1971), page 571], J s f dM is nonneg- 
ative, if and only if 

xf(s) > a.e. (F), 

(2.6) 



mf — J n(x)fp(-,dx)jd\>0, 
(1 A \xf\)p(-,dx)d\ < oo. 



Pick ft G L (5) satisfying (2.4) for t > 0. The following process is i.d.: 
(2.7) X±\l f t dM\ 

with Levy measure in (2.1) given by v = F o jy , where Sx19(s,j;)h 

Z>(a,a;) = x/(.)(«) G IE^ ' 00 ). With Q+ = (0,oo) n Q, the function R in (1.3) 
thus satisfies 



R(u) = 1 A / p(s, 
(2.8) Js V 



SUPt 6( Q+ /t(s) /(1 + CtT) 1 

' d\(s). 



snp tm+ f t (s)+/(l + ctn) 
By (2.5), the process X in (2.7) is H-s.s., if and only if 



(2.9) 



+ J^p^ixa^p^fat^ - 1 -iK(x)a- H J2df at . jp(-,da;)J dA 

does not depend on a > for any choice of n G N, £i, . . . ,i n > and 9i, .. . ,0 n G 
Similarly, X is ff-s.s.s.i., if and only if X(0) = 0, and, with obvious notation, 

JJ^ima~ H (0,f a t + h ~ fh) 

+ Jj^ e ixa-*(6,f a „ h -f h ) _ 1 _ iK ( x)a -B( dj fat+h _ f^p^dx)] d\ 

does not depend on a, h > for any choice of n G N, t±, . . . ,t n > and 
9 u ...,9 n eR. 

Notice that X(0) = 0, if and only if mfo = a.e. (A) and xfo(s) = a.e. 

Example 1. Define a Levy measure p on M with \x\ dp,(x) < oo, 
by /u((— oo, — x)) = r(— x) and /i((x, oo)) = r(x) for x > [/^({0}) = 0], for a 
nonnegative r G L°(R) nL 1 ([— 1, 1]) that is monotone and vanishes at infinity 
on both half-lines. 
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Pick an H > 0. Let M be an i.d. random measure on (0, oo) (equipped 
with the Borel <r-algebra), with Lebesgue control measure, and with 

P (s,B) = Hs-Ms- H B) and m(s) = [ *(*),(., dx). 

JR 

Pick / G L°((0, oo)) satisfying (2.4). Consider the i.d. process X in (2.7), 
where 

(2.10) f t ( s ) = ^ s / t ^ *J>jj> for S >0. 

This process X is H-s.s., since the integral (2.9) evaluates to 

Moreover, we get that X is P-continuous, from the fact that 

E r e *fl(X(t+fc)-X(t))j. 



exp 



s «*(/(»/(t+h))-/(V*)) _ i)^£_^M rfs 



If / = l(o,i]) then X has independent increments, so that P-continuity 
and separability give local boundedness [e.g., Sato (1999), Theorem 11.5]. 

We conclude this section by stating a special case of Rosihski and Samorodnitsky 
[(1993), Theorem 2.1], that is sufficient for our needs, for easy reference. 

Theorem A. Let {X(t)}t£T be an i.d. stochastic process with no Gaus- 
sian component, and with Levy measure v given by (2.1). Assume that the 
parameter space T is countable, and that 

p(sup|Z(t)| <oo) = l. 

// the function 

H{u) = lAv(\y£ R T :supy(t) > u 
is subexponential, then we have 

lim — — -rV\$xrpZ{t) > u \ = 1. 
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2.4. Representation of H-s.s. a-stable processes. Let X be strictly te- 
stable H-s.s., with a E (0, 2). In the case a = 1, assume in addition that X 
is SctiS. 

Let wo 6 L°(5) be positive, and pick a constant /3 £ [—1,1] ((3 = if a = 1). 
Let M be an i.d. random measure (see Section 2.3), with control measure 
A, and with 



p(s, dx) = W (s)((l - /3)l(-oo,0)(£) + (1 + P) 1(0,00)0*0) 



dx 



(2a)\x\ a + 1 ' 



m(s) 



k(x)p(s, dx), 



[k(x) — x)]p(s, dx), 



if a < 1, 
if a = 1, 
if a > 1. 



We say that M is a strictly a-stable random measure. [It is an exercise to 
deduce from (2.2) that M(A) is strictly a-stable (SaS if a = 1), for A G ©o-] 

There exists {/t}t>o ^L> a (5, woA) [which is what (2.4) reduces to here], 
such that X satisfies (2.7), for some (3 € [— 1, 1] (e.g., f3 = ±1 works if a ^ 1). 

Now a process given by (2.7) is strictly a-stable. Denoting a;W = \x\ a sign(x), 
(2.5) shows that X is H-s.s., if and only if /o = a.e. (A), and the following 
integrals do not depend on a > for any choice of n € N, ii, . . . , t n > and 
01, . . . , 6 n € R: 



-if 



i=i 



wo ^A 



(n \ ( Q ) 



Further, X is H-s.s.s.i., if and only if /q = a.e. (A), and the following 
integrals do not depend on a,h > for any choice of n G N, ii, . . . , t n > 
and 6q,...,0 n eR: 



^ 5Z Ojifatj+h — fh, 



wo <iA 



and 



/ f^5>(/< 

J5 V .7=1 



atj+h 



fh 



wo d\. 



Remark. Much is known about the class of .ff-s.s.s.i. a-stable processes, 
that is very rich for a < 2, unlike the Gaussian case. See, for example, 
Samorodnitsky and Taqqu (1994), Surgailis, Rosinski, Mandrekar and Cambanis 
(1998), Burnecki, Rosinski and Weron (1998) and Pipiras and Taqqu (2002a, 
b). 
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For H £ (1/a, 1] with a > 1, and for H = 1/a > 1, it is known that (sep- 
arable) -ff-s.s.s.i. a-stable processes are locally bounded. For other values 
of H and a, local boundedness is not determined by H and a, and there 
exist both locally bounded and unbounded processes. Precise conditions for 
local boundedness are known for a < 1. See Kono and Maejima (1991) and 
Samorodnitsky and Taqqu (1990, 1994). 

3. Overload and the Piterbarg properties. Here we first study the probabil- 
ity for overload P{Y(t) > u}, and then the Piterbarg properties (1.2). 

The next assumptions limit the effect of the left tail of X on the right tail 
of Y: 

, . , P{X(1) < -e(t(«)«-V7)-ff u i-ff/7} 

(3.1) limsup 1 v ; y ' < oo for all e > 0, 

(3.2 lim 1 v ; v \ '- „, i J -L = for some e > 0. 

Here (3.1) is used together with the growth condition indicated in the In- 
troduction 

t(u) 

(3.3) limsup —, < oo. 

In (3.2), v:IR^ (0,oo) is a suitably selected function, with the following 
properties: 

v(u) t(u) 

(3.4) limsup < oo and lim — — u , ztttstt^ = 0- 

In the hypothesis of Theorem 1, (3.1) follows from (3.3) [by (3.11) and 
OR], while (3.2) follows from (3.4) [by (3.20) and R£ PD], provided 

that 

fo^ v Pffll)<-«} 

(3.5) limsup— - r— ■ T <oo. 

u->oo P|SUp s6[0il] X(s) > Uj 

Notice that (3.5) holds if, for example, X(l) is symmetric or nonnegative. 

Theorem 1. Let X be H-s.s. and i.d. with no Gaussian component, 
satisfying Condition X. Consider the process Y , given by (1.1), with c > 
and 7 > H constants. Suppose that the function R, given by (1.3) [with v 
given by (2.1)], belongs to S n PD. Then Y(t) < oo a.s. for each t>0, and 



(3.6) 
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(i) If X is s.L, then for t(u) > 0, 

(3.7) Urn P <7'"» >">=!. 

(ii) If i? E OR and t(u) >0 satisfies (3.1) and (3.3), then we have 

, 00 , n . .P{Y(t(u))>u} _ p{y(t(u))>«} 

(3.8) < hm mf / ) '' - < lim sup , > ^ - < oo. 

(iii) If there exists a function v satisfying (3.4), s-uc/i t/iai 

(3.9) lim = l, 

and such that (3.2) holds, then (3.7) holds. 

Proof. Let ReSf] PD, and denote c = 2 7 /c. For every i > 0, we have, 
for u > 0, 

p{y(t)>n}<p|x(t)<-|| 

+ p{ sup X(s)>^\+p(snp(X(s)-syc)>^\. 

Since X is locally bounded, the first two terms on the right go to zero as 
u — > oo. Furthermore, we can bound from above the third term by 



+ £p su p x 

= P sup X(s)>^—\ 

+ sup X( S )>± ? 

i=0 U G [1/2,1] C J 

^ r x( s ) 2(7- fl )j« 1 - ff /7 

<2VP sup — >^ T7 r- 

" ^ U[o,i]l + caT (l + c)(2Vc) 

Here i? E 5, together with Theorem A, gives 



2™u 

(s)>-r- 



(3.11) 



lim sup r P< sup X(s) > (1 + c)itl 
u^oo it(u) Ue[o,i] J 

< lim sup ^r-^Pi sup — > u } < 1. 

n^oo /t(u) Ue[o,i] 1 + cs 



Here and in future applications of Theorem A, we use that the process under 
consideration is separable and P-continuous. Hence it is enough to consider 
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suprema over any countable dense subset of the parameter space of that 
process [e.g., Samorodnitsky and Taqqu (1994), Exercise 9.3], which we take 
to be the rational numbers in the interior of the parameter space (when that 
parameter space is an interval). 

Our R £ PD has a so-called upper Matuszewska index a < [e.g., Bing- 
ham, Goldie and Teugel (1987), page 71] such that, given —a £ (o, 0) and 
A >0, 

-^-r- < C A~ a for A > Ao and u large enough, 
R(u) 

for some C > 0. Hence the right-hand side of (3.10) is at most 

'2(rf- H )j u 1 - H /"t 



2 y^2R[ - - 



< 4C(1 + c) a (2 V c)° I ^ 2~ a ^~ H ^ J R(u l ~ Hh ) 
\j=o 



for u large enough. This proves the fact that Y(t) < 00 &.s. □ 
Further, by self-similarity and Theorem A [cf. (3.11)], 

p{y(o) >u} = p/sup-^- > 

(3.12) [ s >o 1 + cs 1 

as u — > 00. 

Proof OF (i). For X s.i., y is stationary, and so (3.7) is the same thing 
as (3.6). □ 

Proof of (ii). By (3.3), we have, for some 6 £ (0,1), for all u large 
enough, 

(3.13) 9< mf 1+C '-"" 1/7 ) 7 < sup l + c( a - to -VT r<1 

~ a >t«-V7 1 + CS7 ~ s > tu -i/7 1 + CST 

[where t = t(u)]. Using self-similarity, we therefore obtain, for u large enough, 

U>tu-V7 l + c(s-tU V7)7 



(3.14) 



>P^ sup X U-?»- lh K u^\, 
,>tu-V-» 1 + cs 7 
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Notice that, denoting 

m ( u ) = {tu~ lh )~ H {l + c(Ktu- 1 hy), 
m {u) = {tu~ lh )~ H {l V (Ktu-^y, 

for a constant K > 1, we obtain 

>p( sup X(s)-X(tu-y-r )>ul _ Hh 

> P ( sup J^L >2u 1 - H /A 

- P^tn" 1 / 7 ) > (1 + c(Ktn- 1/7 ) 7 )n 1 -^ 7 } 

(3.15) -P{I(l)>t )1 ( U )ii 1 - |, /' f } 

>p(sup^> 2 ^^;" g/7 i 

- U>Yl + CS 7 J 



Pi sup- 

U>0 1 



+ CST " 1 + C 



Here we have, picking a constant L > 1, 

X(s) 
sup > u 

s >0 1 + CST 

<p( sup >ul +p( sup X ^ >U 

U>L-i! + CS 7 J lo<5<L-! 1 + CS^ 

sup i i \ >u} + Pi sup — ; - > u 

s >i 1 + c{s/Lp J U< s <i 1 + csT 1 + c(s/Lp 

n ( X(s) u ) f X(s) L H u 
< Pi sup ^- > = \ + P<^ sup — > 

It follows from (3.12) and the fact that i? 6 PZ? that, if L is large enough, 
then 

r x{ s ) l h u\ / r z( s ) 

hmsupP sup v ; > — — WP^sup v 7 



1 

> It > < -. 
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Fixing L such that this relation holds, we get immediately 

™f X(s) u 1 If X(s) 1 
P<^sup — > tt > > -P^sup ^- >u} 

for large enough. Therefore, by (3.12), 

(3.16) liminf — 1 „ , p(sup X ^ >u\>l. 

Since i? G PD, and lim sup^^^ r\2 (w) /rji ( u ) < 00, we get from (3.12) 
and (3.16), 

r ■ fp/ £(g) •> 2 ^(")« 1 ' JJ/1 
hm ml P < sup > fj 

u^oo U>Y 1 + csT K H 

(3.17) * {SoT+csi * ^T+c } 

> I lim inf „ f^W-" > ') / a f'W»'- g/T 

- 2 M^°c V # H / ' V 1 + c 

> 2 

for all K large enough [where rji(u)u 1 ~ H /" 1 — > 00 for i = 1,2]. Fixing X > 1 
such that (3.17) holds, we may apply (3.16) and (3.17) on the last row of 
(3.15), to get 

1 T,f X(s)- X{tu- l h) HH 
liminf — ... - Pi sup — ^- - > u l ~ Hh 

J v ■ 1 pf g(g) ^ 2r /2 (^)u 1 - ff ^ 
> - lim ml — , P < sup > 



2 u-oo R(u 1 ~ H h) U>il + ca7" ifH 

1 , p i?(2if- H L^ H ri2(u)u 1 ~ H ^) 
> - hm ml — — F7 — 

- 4 «^oo R{u l ~ H h) 

>0, 

using it! G OR for the last inequality. By (3.14), this gives the lower bound 
in (3.8). 

The corresponding upper bound in (3.8) follows from (3.14). This is so 
because 

limsup \ p(sup X(s) - x( ^" 1/7) > eu^h 

1 X ( S ) . 0..1-IT/7 



< limsup — — : — ttt— P< sup > -u 

~ u^oJ R(u l ~ H h) U>ol + csT 2 

+ hm sup — , \ „\ '- ' . 

t*-oo Riu^h) 
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which is finite, by (3.1) and (3.12), together with self-similarity and R € OR. 

□ 



Proof of (hi). We have 



(3.18) 



1 - 0{{tu^f M ) < inf 

< sup 

< 1 



i + c( s - tu-viyi 

1 + csT 
l + cjs-tu- 1 ^ 

1 + CST 



as u — > oo. Here (3.4) shows that, with obvious notation, 



<o{v{u l - Hh )/u l - Hh ). 
This gives us the following version of (3.14), that for u large enough: 



P{Y(t) > u} < 



1 +CST 



< P<j sup 

s>0 



s>tu~ 1 /-' 

X(s)-X(tu~ l l 



>«i- ff /'r_ et ,( u i-^/7)l. 



(3.19) 

To bound the ratio in (3.7) from below, use self-similarity, (3.4) and (3.12), 
to get 



(3.20) 



'{ sup X(s) >ev{u l - H ^)\ 
< P<^ sup , . '- > 



>0l + CST {l + c){tU- l h) H 
- U>ol + CST 

~ R{Ku l ~ Hh ) 



as u — > oo, 



for any constant if > 1. Hence (3.19), together with (3.9) and (3.12), give 
that 



P{Y(t) >u} 

R{u l - H h) 



lim inf 



> lim inf 



X{s) 
o 1 + csi 



V{X{tu- 1 l~<)>ev{u 1 ~ H h)} 



lim sup ■ 



R{v}~ H h) 
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-limsup ^ , P( sup X(s)>ev(u 1 - H ^)\ 
u^oc R{u l H h) l se[0)tu -l/7) J 

R(Ku l ~ H h) 
> 1 — 2 limsup — — 777— >1 as K — > oo, 

since R G Pi). Of course, this establishes that 

R(u l ~ H h) ~ 

On the other hand, since (3.9) and monotonicity of R give P(u — \v{u)) ~ 
P(u) for any AgR, (3.19) together with (3.2), (3.9) and (3.12), show that 

p{y(t)>4 

Inn sup — : — „, - 

w ^oo R{u l ~ H h) 

< limsup - / „. x p(su P X ^ > u l ~ Hh - 2ev{u l - Hh )\ 

P{-X(l) > e{tu^ 1 h)- H v{u l - H h)} 
+ hm sup — — \ , jj t — 

R{v?-- H h) 

= 1+0. □ 

To establish the Piterbarg property, we make use of the following assump- 
tions: 

P{inf s6[0j i] X(s) < -e(t{u)vT 1 h)- H u 1 - H li} 

hm sup — — — < oo 

(3 21) R{u l - H h) 

for all e > 0, 

P{inf se[oa] X( S ) < -gft^^-V^-^^i-H^)} _ 

(3.22) ™ flCu 1 "^) 

for some e > 0. 

Assumption (3.21) will be used together with the growth condition (3.3), 
while in assumption (3.22), v is a suitably selected function that satisfies 
(3.4). 

In the hypothesis of Theorem 2, (3.21) follows from (3.3) [by (3.11) and 
Re OR], while (3.22) follows from (3.4) [by (3.20) and monotonicity of R], 
when 

(3.23) hm sup — < oo. 

«->oc P{sup se[0jl] X[s) > u\ 

Clearly, (3.23) holds for X symmetric or nonnegative. Otherwise, (3.21) 
and (3.22) could possibly be verified by Theorem A, for — inf se ro ij X(s) 
subexponential, or by Albin [(1998), Theorem 3 and Sections 8 and 9]. 
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Theorem 2. Let X be H-s.s. and i.d. with no Gaussian component, 
satisfying Condition X. Consider the process Y, given by (1.1), together 
with the function R, given by (1.3), where c> and 7 > H are constants. 

(i) Let R€ S n ORH PD. If (3.3) and (3.21) hold, we have, for < 
i(u) < t(u), 

1 ■ r P { SU Pse[0,t(«)] y ( S ) > U ) 

1 < lim ml 7- 

- u^oo P{Y(t(u)) > u} 

P{sup se[0ji(u)] Y(s) > u} 
< Iimsup , < 00. 

~ u^oo P{Y(t(u))>u} 

(ii) Let RgSHPD. Take t(u) and v such that (3.4), (3.9) and (3.22) 
hold. The process Y has the strong Piterbarg property (1.2c). 

PROOF of (i). It is enough to show the upper bound. Using (3.13) and 
(3.14), we get 



sup Y(s) > u I 
Ue[o,t] J 



X(s)-X(r) 1 ffAv 

= P{ sup sup — \L > u x ~ Hh 

U<r<tu-^h s>r l + c{s-rp 

f X(s) 6u l - H ^^ f . . flu 1 "^] 

<p sup— H-> +P sup -X(r)> 

L> 1 + CS> 2 J Lo^r^tw-V-T 2 J 

(3.24) 

for some 9 € (0, 1). Therefore, self-similarity, (3.12), (3.21) and R G (9-R, give 

P{sup se[0ii ] ^(s) > u} 
Iimsup — ttt—, < 00 • 

u-,00 R(u l - H h) 

Now the upper bound desired follows from (3.8) [notice that (3.21) implies 
(3.1)]. 

□ 

Proof of (ii). Using (3.18) together with (3.4), as in the last paragraph 
of the proof of part (hi) of Theorem 1 , we may readily modify the estimate 
(3.24) to obtain 



P<^ sup Y(s) > u \ 
Le[o,t] J 

< p( sup J^lL > U i-Hh _ £v{u i- 

- U>0l + CST V 

+ p{ sup ~X(r) >ev{u l ~ Hh )\. 
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By application of (3.9) together with (3.12) and (3.22), this shows that 



P{sup se[0 1] Y(s) > u} 

limsup — — m— < 1. 



L SUP - 



On the other hand, 

inf Y(s) > u \ 
»e[0,t] J 



X(s)-X(r) , it/ 
inf sup — — t Y- > u l ~ H ^ 

0<r<tu- 1 /i s >r 1 + C(s — Tj* 

>P{ inf (sup . -X(r) + )>u 1 - fl M 

lo^r^tM- 1 /^ \s>r 1 + C(S — r)T / J 

>p( inf sup X ^ - sup X(r)+ >u 1 ~ H/7 l 

l0<r<t«-!/7 s >r 1 + CS^ <r<ti 1 - 1 /7 J 

and, as was established in the proof of Theorem 1, this gives us 

. P{inf se[0)t] Y(s) > u} 
hmmf 7 ' JTI > 1. 

Hence the strong Piterbarg property. □ 

Here are two easy corollaries to Theorems 1 and 2, that make use of (3.23): 

Corollary 1. Let X be H-s.s. and i.d. with no Gaussian component, 
satisfying Condition X. Let the function R belong to S n PL) , and assume 
that (3.23) holds. 

(i) The Piterbarg property holds for 7 < H + (H A 1). 

(ii) If (3.9) holds for v(u) = (1 V uf, for some (3 G (0, 1), then the Piter- 
barg property holds for 7 < H + (if A 1)/(1 — 0). Hence, if (3.9) holds for 
v(u) = (1 Vu)°, /or eac/i /? G (0, 1), i/ien £/ie Piterbarg property holds for any 
7>£T. 

Proof of (i). By part (ii) of Theorem 2, together with an inspection 
of (3.4), it is enough to exhibit a positive function v, with lim M ^ oc t;(u) = 00 
and limsup^^oo v(u)/u < 00, that satisfies (3.9). This is easy: Let bo = 0, 
b\ = 1 and 

f R(x) 11 

6j + i = inf < u > max(6j, 2i) : inf — > 1 > for i > 1. 

^ a;>u i?(x — i) i J 



„l/7-(l-H/7)/(HAl) 
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Since R £ 5 C £ (see Section 2.1), this is an increasing to infinity sequence 
of finite nonnegative numbers, and we may now choose 

v(u) = i if u £E [6j, for i > 1 and i>(it) = 1 for u < 1. □ 

Proof of (ii). Once again, the result follows from part (ii) of Theorem 
2, by means of checking that (3.4) holds for 7 < H + (ff A 1) /(l — (3), when 
u(u) = («Vl) i3 . □ 

In Part II of Example 5, below we see that the Piterbarg property may be 
absent, when, in the notation of part (ii) of Corollary 1, y>H + 1/(1 — (5). 

Corollary 2. Let X be H-s.s. and i.d. with no Gaussian component, 
satisfying Condition X. Let the function R belong to SPiPD, and assume 
that (3.23) holds. 

(i) The strong Piterbarg property holds for t(u) > such that 

_t(u) 

(ii) If (3.9) holds for v(u) = (l\/u)P, for some (3 S (0,1), then the strong 
Piterbarg property holds for t(u) > such that 

*( n ) _ 

(3 - 25) u^c u l/7-(l-«(l-H/ 7 )/(HAl) " °- 

(hi) If (3.9) holds for every v(u) = o(u), then the strong Piterbarg prop- 
erty holds for each t(u) = o(n 1//7 ). 

Corollary 2 is proved in the same way as Corollary 1. 

Notice that the intervals with length t(u) as in (3.25), for which the strong 
Piterbarg property holds, do in fact shrink with u, unless 7 < H + (H A 
l)/(l-/3). 

4. Discussion and examples. Here we present points of view on the re- 
sults of Section 3. Examples are given, where the input process X is repre- 
sented in the form of a stochastic integral, with respect to an i.d. random 
measure. 

We hrst exemplify that the storage process does not have to be finite 
valued, in general. (Obviously, this does not happen under the assumptions 
of our results.) 

Example 2. For standard Brownian {B(t)}t>o motion, and a nonde- 
creasing function / : (0, 00) — ► (0, 00), by the Kolmogorov-Petrowski integral 
test, 



{B(t) < V2i/(i) for t > T, for some T = T(u>) < 00} 
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is a zero-one event, or in other words, 

P{B(t) < V2tf(t) ultimately as t -> 00} = or 1, 
with the probability being 1 if and only if 

/ exp{— fit) } dt < 00 for some to > 

Jto t 

[e.g., Bingham (1986), page 436]. From this we get that 

/ \ (B(t) 2 , , 3 In In In t , 1 

(4.1) P j < mini H ultimately as i -> 00 j- = 0. 

Consider the following H-s.s. process X: 

(4.2) X(t)=t^exp|exp J for t > 0. 

From (4.1) it follows readily that 

X(t) 

limsup = 00 w.p.l for 76K. 

Hence, the storage process Y(t) in (1.1), with input X given by (4.2), is not 

finite valued for any t > 0. [Incidentally, using the fact that {tB(l/t)} t >o = {B(t)} t >o 
together with (4.1), it can be seen that the process X is not bounded at zero.] 

4.1. I.d. H-s.s. processes. In this section, X denotes the i.d. H-s.s. pro- 
cess given in Example 1, which is assumed to satisfy Condition X. Notice 
that, by (2.6), X is nonnegative, if / is nonnegative (nonpositive) and r is 
zero on (— 00, 0] ([0, 00)). 

Denoting, for s > 0, 

<t+(s) 1 = s tt sup \ K ' ' and <r_(s) 1 = s H sup 



1 + an teQ + 1 + cm 

we have, by (2.8), for u large enough, 

roo fjg 
R(u) =H fi(R\ [-a-(s)u, cr + (s)u]) — 
(4.3) Jo 

= H (r(o+(s)u) +r(-a-(s)u)) 



ds 
s 

ds 

' " " s ' 

Example 3. Let r be symmetric, with r G RV{— p) for some p > 0. By 
(2.2), it follows that X is symmetric. Assume that for some e > 0, 

(4.4) / (a + (s)- {p+£) + a_(s)- {p+£) ) — < 00. 

Jo s 
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Then by (4.3), we have 

r°o (To 

(4.5) R(u)~H (a+(s)- p + a-(sy p ) — r{u) asu^oo. 
Jo s 

Thus part (ii) of Corollary 1 shows that the Piterbarg property holds for 
7 > H, and part (hi) of Corollary 2 gives the strong Piterbarg property for 
t(u) = o( U V7). 



In Example 3, symmetry gives us (3.23), for free. Without symmetry, we 
may still establish the Piterbarg properties, by direct verification of (3.23) 
[or (3.22)]. 



Example 4. Take / nonnegative, and not identically zero. Assume that 



r G RV(—p) and r(— •) G RV(— p) with lim sup — -— < oo, 

«->oo r(u) 



for some constant p > 0. Under the condition (4.4), we see that (4.5) holds 
with er_(s) = 0. Further, as in (4.5) by Theorem A, the limit in (3.23) is 

(7°V' sup f( r y—/n s H P sup f(r)P—) 
\Jo re(s,oo)nQ s i Jo re(s,oo)nQ s J 

r{—u) 
x limsup — — — — < oo. 

u^oo r(u) 

Hence the Piterbarg properties hold in the same way as in Example 3. 



Example 5 — Part I. Pick constants A > and a G (0, 1), and consider 

(4.6) r(x) = g{x)e~ Axa for x > and r(x) = for x < 0, 

where g G RV(p) with p G R. Take / nonnegative, so that X is nonnegative. 

Let <7+(s) take its minimal value at a unique s > 0, where cj+ is two 
times continuously differentiable at s, with o+(s) > 0. By Taylor expansion 
in (4.3), we have 

R(u) = Hu-^ r g(a+(g + S / U t )u) e- A ^ s +^ a/2 ^ ds 



,s 



sy / ^acr^(s)/o-+(s) 1 - Q 
(4.7) 



5 (n)-u- Q / 2 ^exp|-^^ + (s) + i<(s)s 2 /^ Vjds 
g{u)u a/z e 14 asu^oo. 
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Hence, part (ii) of Corollaries 1 and 2 applies, for (3 < 1 — a, to give the 
Piterbarg property for 7 < H + (H A l)/a, and the strong Piterbarg property 
for t > with 

urn — n 7- — „, w7777-r; = for some a > a. 

Example 5 — Part II. Here we continue the study of the case when r 
is given by (4.6), in the particular case when / = l(o,i] (so that X has inde- 
pendent incre- ments). We show that, in this case, the Piterbarg property is 
absent, if 7 > H + 1/a. 

By Theorem 1 and (3.24), the Piterbarg property is absent when 

/ \ , 1 f X(s)-X(r) I 
(4.8) hmsup — r-rP-{ sup sup — — L >u\>l. 

u^oo R(u) Lo<r<t/uV(7-H) s>r 1 + c(s - rp J 

Theorem A does not apply here, since suprema are taken over regions that 
depend on it. However, the arguments for that theorem in Rosihski and 
Samorodnitsky (1993) produce an asymptotic lower bound, for the proba- 
bility in (4.8), which implies the following sufficient condition for (4.8): 

lim sup 



R(u) 

x.LlMnQ; sup sup y(*)-y(r) >u 

{ rg(0,t/-u 1 /(7-H)) n Q se ( r . iOO )nQ 1 + c(s — rp 

> 1. 

(4.9) 

Denoting the numerator in (4.9) by Rt(u), we have, by the inequalities 
7 > H + 1/a and (1 — x) H < 1 — (H A l)x for x E [0, 1], together with (4.7) 
[cf. (4.3)], 

H 1 (o,s](y) - l(0,r](y) ^ \\dy 



R t (u) = H fi[\x>0:y n x sup sup y ^ ' > u 

V I o<r<t/«V(7-») s>r l + c{s-tp ) ) y 

H / r(y u) — 

y 

. roo j i + c ( y _ t / u m>r-H)p ^ dy 



> H 



> H 



t / u i/(7-ff) V y H J y 



oc 



r ( ffto _ t/B v<^ ) (»z^)-.)4 

r (a(y - t/u l ^~^)(l -(HA 1)^) u) ^ 
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s 

X 



~ i?(u) exp|^4(j(s) a (if A asu^oo. 

This gives (4.9). [If nervous about this calculation, shrink the domain of 
integration from \t/u l ^" 1 ~ H \ oo) to s±K/u a / 2 , and send K j oo at the 
end.] 

4.2. a-stable processes. First we consider a storage process Y, with an 
a-stable H-s.s. input process X. 

Example 6. Let X be a strictly a-stable H-s.s. process, that satisfies 
Condition X, and is given by (2.7), where M is a strictly a-stable random 
measure (SaS if a = 1). By calculations similar to those in Examples 3 and 
4, we have 

flW . 1A „-«/(i^ 8up M£)_£ 

V 2 teQ + (1 + cf) a 

+ — ^ + JTTcW) Wo{s)dX{s) - 

We assume that the above integral is nonzero, so that R is not identically 
zero. 

Provided that X satisfies (3.23), Corollary 1 now gives the Piterbarg 
property for any 7 > H, for the storage process Y, while Corollary 2 gives 
the strong Piterbarg property for t(u) = o(n 1 / 7 ). However, by Samorodnitsky 
(1988) [see also Samorodnitsky and Taqqu (1994), Theorem 10.5.1], the limit 
in (3.23) is 



/(i±£ sup (/■-)« + !_£ sup (f+r)w Q dx) 
JS\ * te(o,i)nQ L te(o,i)nQ / / 

sup u r r+ l -±l sup (f+AwodxY 1 . 

JS\ * tefo.ilno z teco.imo / / 



te(o,i)nQ z te(o,i)nQ 

This ratio is finite, by the local boundedness of X and the assumption that 
R is nonzero. 

Next we consider the Piterbarg properties (1.2), in the case when the 
process Y itself is a-stable. Now Y is no longer a storage process, and the 
example simply is to illustrate how unusual the Piterbarg property is for 
"usual processes." 
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Example 7. For Y = {Y(s)} s >o an a-stable process, a 6 (0, 2], we may 
write Y = Y\ — Y2 + fj,, where \i : [0, 00) — > R is a suitable function, while 
{Y"i(s)} s >o and {l2(s)} s >o are independent a-stable processes, such that (for 

3 = M 

and s > 0) 

exp{-a?. (s) |er(l -itanf™) sign(0)) }, if a ^ 1, 



E{e^ (s) } 



exp j-ay j(s) |0| + ^ sign(0) ln(|0|)H, if a = 1, 



for e R. 

Here cry^( s ) is the scale parameter of the a-stable random variable Yj(s). In 
the Gaussian case a = 2, we may take Y\ = 0. We assume that Y satisfies 
Condition X, from which it follows that Y\ } Y% and fx can be taken to satisfy 
Condition X. 

We are going to investigate when the Piterbarg property (1.2a) holds. 

d 

Case 1. If Y x / 0, then (1.2a) holds if and only if the f.d.d.'s of Y\ 
coincide with those of a single a-stable random variable [Samorodnitsky 
and Taqqu (1994), Theorem 10.5.1]. 

Case 2. If Y x = and a > 1, then (1.2a) holds if and only if the f.d.d.'s 
of Y<i coincide with those of a single a-stable random variable, and /i is 
constant. 

To see this, notice that ay 2 and fx must be constants on [0,t] [Samorod- 
nitsky and Taqqu (1994), equation 1.2.11]. Given these properties, we have 



(4.10) 



sup Y(s) > u > 

iGfO.tl J 



se[o,t] 

> P{{Y(r) >u}U {Y(s) > u}} 

> 2P{Y(r) > u} - P{i(y(r) + Y(s)) > u} 

for r, s G [0,t]. If a > 1, then the second probability on the right-hand side is 
o(P{Y(r) > u}), unless &Y(r)+Y(s) = °Y(r) + °Y(s)- By Minkowski's inequal- 
ity, this happens if and only if Y(r) = Y(s) a.s. If a = 1, then the spectral 
measure T of (Y(r), Y(s)) [Samorodnitsky and Taqqu (1994), Section 2.3] is 
supported on = {(si, S2) £ K 2 : Si + s 2 . = 1, si, S2 < 0}, and [Samorodnit- 
sky and Taqqu (1994), Example 2.3.4] 

Y(r) + Y(s) 



Y{r)--{( S -^ln S -l±^dT(s)- [ Sl ln\ Sl \dT(s 
n \JS~ 2 2 Js~ 
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By a convexity argument, unless Y{r) = Y{s) a.s., the term f (•••) on the 
right-hand side is strictly positive, so that the second probability on the 
right-hand side in (4.10) is o(P{Y(r) > u}) [Samorodnitsky and Taqqu (1994), 
equation 1.2.12]. 

See Talagrand (1988) and Albin (1999) for more information related to 
Case 2. 

Case 3. If Y\ = and a < 1, then (1.2a) holds (in the sense of 0/0 = 1) 
since I2 is nonnegative. 

Special case. If Y is SaS, then (1.2a) holds if and only if the f.d.d.'s 
of Y coincide with those of a single a-stable random variable. 

Turning to (1.2b), with t = t(u) — > 00 as u — > 00, the above character- 
izations remain valid [with (1.2b) replacing (1.2a)], if appropriate global 
boundedness properties are imposed on Y2 and /i in Case 1, and on fi in 
Case 3. 
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